A relativistic treatment of pion wave functions in the annihilation pp — > n n 
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Quark model intrinsic wave functions of highly energetic pions in the reaction pp — > tt~-k + are 
subjected to a relativistic treatment. The annihilation is described in a constituent quark model 
with A2 and R2 flavor-flux topology and the annihilated quark-antiquark pairs are in 3 Po and 3 Si 
states. We study the effects of pure Lorentz transformations on the antiquark and quark spatial 
wave functions and their respective spinors in the pion. The modified quark geometry of the pion 
has considerable impact on the angular dependence of the annihilation mechanisms. 
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I. INTRODUCTION 

The LEAR experiments Q on pp — > 7r~n + and pp — > 
K~K + , which yielded a large and very accurate data set 
for differential cross sections da/dQ and analyzing powers 
A 0n from 360 to 1550 MeV/c, have until very recently 
resisted a satisfying comparison with theoretical models. 
In particular, the observable A 0n exhibits a characteris- 
tic double-dip structure which is hard to reproduce. An- 
other prominent feature of the Aq u data is a shift of the 
asymmetry from predominantly negative values at lower 
momenta toward positive values at higher momenta. As 
already mentioned in , this indicates the presence 
of several partial waves while most model calculations 
[J B l e& d to amplitudes dominated by total 
angular momentum J = and J = 1. The reason for this 
feature of all models is a rather short range annihilation 
mechanism,either driven by baryonic exchange in the t 
channel \^ |fj, |fj, |9j or by overlap of quark and antiquark 
wave functions of the proton and antiproton as in [21 la| • 

One way to improve models would hence be to in- 
crease the annihilation range. In Ref. [2|, this prob- 
lem is addressed within the framework of a constituent 
quark model. For a summary of the pp annihilation in 
a quark model see the review by Dover et al. |lfjj |. In 
paper 0, the radii of the proton, antiproton and pion 
were readjusted so they coincide with the ones obtained 
from measurements of the respective electric form fac- 
tors [TTL IT2| rather than with the considerably smaller 
constituent qqq and qq quark core radii. The latter radii, 
which ignore the hadronic qq cloud, were used in earlier 
attempts to reproduce the LEAR observables 0, ■ In 
the model calculations, this increase in radius is directly 
related to the Gaussian description of the intrinsic quark 
wave functions. In the overlap integral over the quark 
coordinates, the wider Gaussian wave functions provide 
in turn a larger overlap of the quarks and antiquarks and 
therefore a larger range of the annihilation mechanism. 
In fact, the radii in Ref. were obtained in a fit to 
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da/dfl and Aq u data mentioned above and are about 7% 
larger than the charge radii derived from the electric form 
factors. This effect, along with a fine-tuned final-state in- 
teraction of the 7r~7r + pair, improves the reproduction of 
the observables da/d£l and A 0n , and in particular of the 
double-dip structure of the latter one, considerably. 

As was pointed out in Ref. , there is another reason 
to believe that the geometry of the final pions is rele- 
vant to higher partial wave (J > 1) contributions to the 
total amplitude. At the center-of-mass (cm.) energies 
y/s considered in the LEAR experiment pp — ► 7r~7r + , fi- 
nal pions are produced with kinetic energies much larger 
than the pion rest mass. For example, for an antipro- 
ton beam with pj ab = 800 MeV/c the total cm. energy 
is E cnl = y/s ~ 2020 MeV. This translates into a rela- 
tivistic factor 7 = i? cm /2m w = 7.2. In the quark model, 
the pions are described by Gaussian spheres in the rest 
frame of each pion. Yet, the transition amplitudes and 
implicitly the observables of the pp — > tt~it + reaction are 
calculated in the cm. frame. In the cm. frame the 
outgoing pions, due to their large kinetic energy, will be 
shaped like pancakes. This will alter the angular depen- 
dence of the annihilation mechanism. 

This modification of the pions, due to Lorentz trans- 
formed Gaussian wave functions, was investigated previ- 
ously by Maruyama et al. in the annihilation of 
atoms at rest into two mesons. The authors of Ref. [L 
found that the two-pion production is enhanced upon 
including Lorentz contraction effects. However, only 
branching ratios were discussed and the impact of these 
Lorentz contractions on the angular dependence of the 
annihilation pp — > tt~tt + at positive energies remained 
unknown. 

In this paper, we will show that the pion deformation 
affects the overlap of intrinsic pion, proton, and antipro- 
ton wave functions. It was in this spirit that in Ref. 
the width of the Gaussian wave function was increased 
not merely to reproduce the measured charge radii, but 
also to mock up an altered overlap of the quarks and 
antiquarks implied by relativistic considerations. In the 
following, we will carry out the Lorentz transformation 
which is a twofold task: the pions are qq pairs described 
by the usual spinors of the free Dirac equation times a ra- 
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dial Gaussian function to account for their confinement 
within the pion. Hence, the Lorentz transformation is 
effected on both spin and Gaussian components of the 
wave functions. Using the relevant Feynman diagrams, 
we analyze the relativistic effects on the amplitudes and 
compare them with previous results Q- 



II. LORENTZ TRANSFORMED PION 
INTRINSIC WAVE FUNCTIONS 

The reaction pp — > tt~tt + is described in the cm. 
frame. Therefore, given the large kinetic energy of the 
final-state pions, we must Lorentz transform each pion 
intrinsic wave functions from the pion rest frame to the 
cm. frame. As mentioned before, this will affect the 
Gaussian radial part of the wave function which will be 
distorted along the boost direction into an ellipsoid. We 
begin with an r-space realization of the pion wave func- 
tion and first concentrate on its spatial component. The 
proton and antiproton are described by similar Gaussian 
functions but given their much smaller kinetic energy and 
larger mass their intrinsic distortions are ignored. The 
antiquark and quark coordinates in each pion transform 
according to 
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(i) 



We split the quark and antiquark coordinates into com- 
ponents parallel and perpendicular to the boost direc- 
tion (3 = v/c, which is related to relativistic boost factor 
7 = (f — /3 2 )~2 = E cm /2m n . Hence, the parallel compo- 
nent of the separation — rj between an antiquark-quark 
pair in the pion rest frame is related to the one in the 
cm. frame with coordinates y[ and by 



3>\\' 



(2) 



Here, we have used the equal time condition t[ = t[j in 
the cm. frame. These transformations can be applied 
straightforwardly to the spatial part of the pion wave 
functions. Inserting the Lorentz contraction in Eq. J2J, 
the pion wave function suppressing spin, isospin, and 
color dependence, 



tp^T^Tj) = A^expj - - ^(rv-R^) 2 }, 

a=i,j 



(3) 



becomes 



N, 



ex p{"f E [(r Q -ar)!+7 2 (r tt -R^]} (4) 



a=i,j 



where f3 is in this case the size parameter and = 
\{vi +Tj) is the pion coordinate. The new normalization 
factor N n = J^j comes from the condition that Eq. Q 
be normalized to unity in the cm. In Eq. ijljl. the spatial 
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FIG. 1: a) Rearrangement diagram R2 and b) Annihilation 
diagram A2. The numbers with bars denote antiquarks, those 
without bars the quarks while the dashed lines represent ei- 



ther the exchange of the effective 
3 Si state. 



'vacuum" Po or "gluon" 



distortion of the pions in the cm. of the reaction pp — > 
7r~7r + alluded to in the introduction is shown explicitly. 
In p-space the wave function is 

x ex P {" Yd E [( p « " i p -l 2 +i(P« - } 
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The momenta Pi and Pj are the quark and antiquark mo- 
menta, while P^ = pj + pj is the pion momentum which 
is parallel to the boost direction. With these Lorentz 
transformed wave functions, we are now in the position 
to evaluate the 3 Si and 3 Pq amplitudes taking into ac- 
count the boosted pions. We will consider two types of 
commonly used diagrams, the rearrangement R2 and an- 
nihilation A2 diagrams as depicted in Fig. ^ 



III. EFFECTS OF DEFORMED PIONS ON THE 
ANNIHILATION AMPLITUDES 



With the Lorentz transformed pion intrinsic wave 
functions.we can recalculate the T matrix elements of 
Refs. 0, Q for the annihilation reaction pp — > 7r~7r+. 
Here, we will do this first for the R2 diagrams (which we 
used exclusively in a previous paper 0). The transition 
operators T( 3 Pq) and T( 3 Si) are obtained from integrat- 
ing out the quark and antiquark momenta. In the fol- 
lowing, P p , Pp, P.„-- and are the p, p, ir~ , and ir + 
momenta, respectively. For the R2 diagram of Fig. [TJi, 
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the transition operators for 3 Po as well as 3 Si are obtained from the integral: 
I 



T R 2( 3 Po, 3 Si) = N J dpi dp 2 dp 3 dp 4 dp 5 dp 6 d Pl dp' 2 dp 4 dp;, S(p[ + p 5 - P.- ) S(p' 2 + p' 4 -P n+ )x 

X S(Pl + P2 + P3 - Pp) £(P4 + P5 + P6 ~ Pp) $(PG + P3 + Pi - Pi) S(p' 2 - p 2 ) S(p' 4 - p 4 ) S(p' 5 - p 5 ) X 

X eXP {-Jf3 



x exp 





(Pi-P5r + (p 2 -P4r+ -^-1 (p' 2 -p 4 )-p. + (pi-p' 5 )-p. 



^ 2 ( 3 Po, 3 ^i), 



(6) 



r 



where we have absorbed the relativistic factor 7, the nor- 
malization factor Nn , and the size parameter j3 in Eq. 1|5|1 
(and similarly N p , Np, and the size parameter a in the 
pp wave function) into an overall normalization Af. 

The two mechanisms differ only in the operators 
Vr2( 3 Po) an d Vr 2 ( 3 Si) stemming from the respective an- 
nihilation mechanisms: 



V R2 ( 3 P ) = er 3er (p6-p 3 ) 



(7) 



Vi? 2 ( 3 5i) = 2<t 36 -Pi + i(<T 3 6Xo- ia )-(pi -Pi')- (8) 

The exponential parts in Eq. 10 originate in the p- 
space component of the boosted pion wave functions 
of Eq. (JSJ, and the baryonic (non-boosted) wave func- 
tions, respectively. The spin-momentum dependency in 
Eqs. (0 and (JHJ) comes from the computation of the R2 
diagrams using the Dirac spinors. The indices on the 
Pauli matrices denote the spinor subspace in which they 
act. For instance, the antiquark-quark pair q^q^ is anni- 
hilated and momentum p 3 + p6 is transferred from this 
vertex to quark qi . 

After Fourier transform into r-space, Eq. ^ becomes 
for the 3 Pq mechanism 

f R2 { 3 P a ) = iN[A v a-R' + BycrR + CV(er-R') i?cos 9] 
x cxp{ AR' 2 + BR 2 + CRR' cos 9 + BR 2 cos 2 9}. (9) 

The various coefficients will be given below. Since the 
3 Si mechanism in Eq. © contains two terms, we split 
their contributions conveniently into a longitudinal and 
a transversal part. It will be seen that they give rise to 
different selection rules. The transition operator is thus 

f R2 ( 3 S 1 ) = N[A L icr-R' + B L io--R + C L icr-~R' R cos 9 
+ A t ct-R' + B t ct-R + C t (ct-R')Rcos9] 
x exp{ AR' 2 + BR 2 + CRR' cos 9 + DR 2 cos 2 9}. (10) 



The relative tt~ tt + and antiproton-proton coordinates are 
R' = Rjr-— Rtt+ and R = Rp — R p , respectively. The 
angle 9 is between the relative ir~ir + coordinate R' and 
antiproton-proton coordinate R in the cm. system. If 



compared with the results of Ref. one observes the 
appearance of additional terms in both transition oper- 
ators, namely the D term in the exponent and the Cy, 
Ct, and Ct terms. These new terms introduce mani- 
fest angular dependence. More precisely, the relativistic 
pion wave functions lead to additional angular depen- 
dence in the form of quadratic cosine terms in the ex- 
ponentials and linear cosine terms in the coefficients of 
the spin operators. Note that aside of the explicit non- 
locality, the transition operators T R2 ( 3 P ) and T R2 ( 3 Si) 
are now energy dependent via the Lorentz boost factor 
7. This 7-dependence is non-trivial and affects the an- 
gular dependence of the transition operators in the sense 
that the various terms mentioned above behave differ- 
ently for increasing values of 7. Also, the magnitude of 
these terms differs immensely depending on whether 7 is 
small (7 ~ 1 — 2) or large (7 — 5 — 8). For large 7 the D 
term dominates the exponential part and, for example, 
Cy ^> Ay, By and Ct ^> At,Bt, however, Ct, At, and 
Bt are of the same order of magnitude. 

The above results for the transition operators are valid 
for the diagram in Fig. ^1. Summing over all possible 
permutations of Fig. where for example p 3 + P6 is 
transferred to quark q 2 or one of the antiquarks (74,^5, 
yields the same symmetry properties as in Ref. Q , as ex- 
pected. In particular, the vacuum coefficients Ay, By, 
and Cy add up such that Tr 2 ( 3 Po) contributes only to 
l 7rn = 0,2,4,... and lp p = 1,3,5,... and hence acts in 
pp states with J n = + ,2 + ,4 + ,... waves. The longi- 
tudinal part of the transition amplitude Tr 2 ( Si), with 
the coefficients At, Bt, and Ct, bears the same symme- 
try properties and therefore acts in the same waves as 
the "vacuum" mechanism. The transversal component 
with At, Bt, and Ct, on the other hand, contributes to 
l 7rn = 1,3,5,... and lp p = 0,2,4,... and acts therefore in 
pp states with J 7 * — 1~, 3~, 5 _ , ... waves. 

The complete form of the R2 transition operators for 
the vacuum 3 Po amplitude when summing over all per- 
mutations is 
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T%-('Po) = 

= iN\A v a- •R'sinh(CR-R') + Sycr-Rcosh(CR-R') 

+ C v (<x-R')i?cos(9cosh(C'R-R' 

x exp{AR' 2 + BB 2 + DR 2 cos 2 8}. (11) 

The total 3 <Si amplitude for the longitudinal component 
is given by 

fg-( 3 Sf) = 

= iN \A L er ■ R sinh(C* R R) + S L cr-Rcosh(CR-R') 
+ Cz,(o--ROflcos0cosh(C'R-R / 
x exp{AR 2 + BR 2 + L>R 2 cos 2 0} 

and 



(12) 



M Atct-R cosh(CR-R') + B T er-Rsinh(C*R-R) 
+ C T (er-R').Rcos6>smh(C*R-R)] x 
x exp{AR' 2 + BR 2 + DR 2 cos 2 0} 



(13) 



for the transversal component. 

The explicit expression for the coefficients of Eqs. (|llfl . 
|Q, and (E3 is: 



A = 

B = 

C = 

D = 

A v = 
B v = 
C v = 



a(5a + 4/37 2 ) 
' 2 (4a + 3/3 7 2 ) 

3(7a 2 + 18a/? 7 2 + 9/3 2 7 4 ) 
8 (4a + 3/3 7 2 ) 

3a(a + /37 2 ) 
"2(4a + 3/?7 2 ) 

9/3( 7 2 - 1 



D 



(5a + 4/?)(5a + 4/3 7 2 ) 



a(a + /?7 2 ) 
4a + 3/3 7 2 

3(a + /j7 2 )(5a + 3/3 7 2 ) 

2(4a + 3/3 7 2 ) 
3/3( 7 2 -l) 



C v 
a 2 



(5a + 4/3)(5a + 4/3 7 2 ) 



(14a) 
(14b) 
(14c) 
(14d) 

(15a) 
(15b) 
(15c) 



A L = 
B L = 

C L = 



-A v 

9(a + /3 7 2 ) 2 
2 (4a + 3/3 7 2 ) 
3/3( 7 2 - 1) 



-Cr 



(5a + 4/3)(5a + 4/37 2 ) 



(16a) 
(16b) 

(16c) 



A T 
B T 

Ct 



-2A V 
3a(a -\ 



h 2 ) 



Aa + 3/3 7 2 
3/3a 2 ( 7 2 



Ct 
-1) 



(5a + 4/?)(5a + 4/3 7 2 y 



(17a) 
(17b) 

(17c) 



In the non-relativistic limit 7 — > 1, the results of 
Ref. are recovered as expected. In particular, the 
coefficients D, Cv, Cl, and Ct vanish in this limit. 



Secondly, we consider the A2 diagrams, where the 
quarks and antiquarks are integrated out similarly to 
Eq. © but where the different flavor-flux topology im- 
plies a different momentum transfer. Effectively, this re- 
sults into one delta function less than for the R2 dia- 
gram. This will affect the angular dependence and selec- 
tion rules as will be seen shortly. 



We take into account two types of A2 diagrams — one 
in which both q^ and q^q2 pairs in Fig. \T]p annihilate 
into a 3 Si state, while in the other diagram the qgqs pair 
annihilates into the "vacuum" 3 Pq state and the q§q2 pair 
annihilates into the "gluon" 3 5i state followed by the 
creation of a q^ qi> pair. We do not take into account the 
double 3 Pq annihilation since the resulting operator is of 
the order C[(p/m) 3 ]. The integrals for the two cases are 
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t A2 (^ , = n J d Pl d P2 d P3 d P4 dp 5 d P6 d P ; dp' 2 dp^ dp^ 5 (p; + p^ - P ff - ) *(p' 2 + p^ - P ff + ) X 

X (J(pi + p 2 + p 3 - Pp) <5(p 4 + P5 + P6 - Pp) ^(P6 + P3 + Pi Pi) <KP2 + P 5 ~ P2 - P5) <5(P 4 - P4,) X 



xex Pi -- 



x exp 





(p'i - Ps) 2 + (Pa - P4) 2 + -o - 1 (Pa - P4)-P- 



i=4 



(Pi-Pe)-Pir 



(18) 



r 



The two cases differ only in the operators Vai according 
to the type of annihilation mechanisms used: 



Vr 2 ( 3 P ) 



X <T25'CT2'5' 



V A 2(^) = V R2 ( 3 S 1 ) x* 25 -v 2 , 5 , 



(19) 



(20) 



The notation is the same as for the R2 diagrams. Fol- 
lowing Fig. momentum is transferred from the anni- 
hilated qeq3 pair to quark q\ (but can also be transferred 
to qi) and from the q§q 2 pair to the created pair q~yq 2 '- 
In r-space the A2 transition operators for the mixed 3 Po 
and 3 5i vertices are 



T A2 (^) = Af[(T-V R > + (A v + B v )ia-H]x 
x 5(3R/2 - R') exp{(B + D)R 2 } 

and for two 3 «Si annihilation vertices 



(21) 



TA2(l S 



= J\f[<r-V R i + (A L +B L )ia--K + 

+ (A T + B T )a-R] x 

x <5(3R/2-R') exp{(B + D)R 2 }. (22) 



For convenience the relativistic terms By, Bt, Bt in 
the spin coefficients and D in the exponent are written 
separately from the corresponding non-relativistic terms 
Ay, Al, At, and B. All are listed explicitly below. 

It is a striking feature of the A2 amplitudes that they 
depend solely on the relative proton-antiproton R vector 
but not on the relative orientation of R and R'. This 
result is an implicit consequence of Fourier transforming 
the A2 amplitudes to r-space which yields a delta func- 
tion <5(3R/2 — R'). One does obtain additional terms 
due to the relativistic corrections in the pion wave func- 
tions, yet these are simpler than for the R2 diagrams and 
they do not introduce extra angular dependence. Again 
certain terms, namely D, By, Bl, and Bt, are energy 
dependent via the boost factor 7. In particular, the ex- 
ponential factor D is very large (for 7 ~ 7 it is about 
fifteen times larger than B) and acts like a short-range 
angle-independent cut-off. 

Since the A2 diagrams were not discussed in previ- 
ous work 0, we take the opportunity to analyze their 



symmetry properties. Taking into account spin-flavor 
and color matrix elements, the amplitudes in Eqs. (|21() 
and H22|) are identical for the permutation of the A2 dia- 
gram (Fig. (Ip) in which a 3 Pq or 3 5i state resulting from 
the qeq^ annihilation is exchanged with q~4 rather than 
with q\. This means that the complete form of T^' re- 
mains as in Eqs. (|21|l and ll'L'li and that Ta 2 (^) and 

Ta 2 {^ l ) are odd in R and R'. Because of the delta 
function <5(3R/2 — R') all A2 operators, regardless of the 
annihilation mechanisms employed, act in (l n7r ^= lp p ) 
,r = O+jl-,2+,3-,4+,5-,... waves. 

We end this section with the expressions for the coef- 
ficients of the transition amplitudes (|21(l and Q22[): 



B = - 



3(28a 2 + 36a/3 + 9/3 2 ) 
8(4a + 3/3) 



D = -9/3(7 2 - 1) -h 



4a 2 



(4a + 3/3)(4a + 3/3 7 2 ) 



A x 

B v = 

A L = 
B L = 

A t = 

Bt = — 



3(12a 2 + 14a/3 + 3/3 2 ) 

2(4a + 3/3) 
3/3( 7 2 - 1) 



1 



4a 2 



(4a + 3/3)(4a + 3/3 7 2 



3(2a + /3) 
2 

3/3( 7 2 - 1) 



6a(a + /3) 
4a + 3/3 
6a 2 /3( 7 2 



1) 



(23a) 
(23b) 

(24a) 
(24b) 

(25a) 
(25b) 

(26a) 
(26b) 



(4a + 3/?)(4a + 3/3 7 2 ) ' 

Clearly, D, By , Bl, and Bt vanish in the non-relativistic 
limit 7 — ► 1. 



IV. GEOMETRIC INTERPRETATION 

It is important to understand the geometric implica- 
tions of deformed pion wave functions as required by rel- 
ativity and their relation to the pp annihilation range. In 
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Sec. 0] we recall that in a previous (non-relativistic) at- 
tempt to explain the LEAR data on pp — ► tt~tt + an 
increase of proton, antiproton, and pion radii effectively 
augments the range of the transition operators Tr2 and 
Ta2- It is therefore surprising that one can obtain any 
improvement from the Lorentz contraction of pions jlSj . 
After all, we "shrink" the intrinsic pion wave function 
in the boost direction, which amounts to less overlap of 
the quark and antiquark wave functions. The graphic 
description in Fig. [3 visualizes this. The two intrinsic 
pion wave functions are plotted at a fixed distance |R'| 
from each other and the centers of the (also Gaussian) 
proton and antiproton intrinsic wave functions are sep- 
arated by a distance |R|. One sees that when the pion 
wave functions are Lorentz contracted in the cm., the 
overlap of pions, proton, and antiproton (and of a transi- 
tion operator depending on the annihilation mechanism, 
which we omit here for simplicity) decreases consider- 
ably and so does their contribution to the annihilation 
amplitude. Since the total probability is conserved, the 
Gaussian wave functions of the pions are peaked higher 
in the lower graph. 

As discussed in Sec. IIIII the R2 transition operators 
acquire new terms which introduce additional angular 
dependence. The effect of the D term in Eqs. 10 and 
(fTOft . for example, is strongly angle dependent — for R 
parallel or anti-parallel to R', the R-range of Trz, is 
drastically reduced due to large D values, whereas for 




FIG. 2: The upper graph shows a two-dimensional projec- 
tion of the intrinsic pion wave functions (solid lines) given by 
Eq. @ and (anti) proton wave functions (dashed lines). The 
distance between the peaks of the solid-line Gaussians repre- 
sents the relative pion coordinate R' = R,r — R,r+ whereas 
the distance between the peaks of the dashed-line Gaussians 



is the relative pp coordinate R 



The distances 



are the same in both graphs, however, in the lower graph the 
boosted intrinsic pion wave function is given by Eq. Q with 
7 = 7. 





Relativistic 



Non-relativistic 



FIG. 3: Sketch of the overlap in the reaction pp tt~tv + 
in the cm. system at two different angles 9 a . m . between the 
incoming pp and the outgoing nn pairs. In the left column, 
the pions are Lorentz contracted, while in the right column 
they are described by spherical wave functions. One sees that 
the angular dependence of the overlaps at same distances R' 
and R strongly differs in relativistic and non-relativistic case. 



R 1 R' this cut-off effect is absent. This angular de- 
pendence is schematically depicted in Fig. |3[ where R' 
and R have fixed magnitudes. For R _L R' in the top 
panels one sees that even though the pions are Lorentz 
contracted, each pion still overlaps with the pp pair as in 
the non-relativistic case. This overlap vanishes at small 
angles (see lower panels in Fig. Ell an d at 8 — (8 = ir) 
when R' is (anti)parallel to R. This illustrates that the 
angle dependence of the overlap in the relativistic and 
non-relativistic case are quite different. That the least 
overlap occurs for R |j R' is expressed in the transition 
operators Tri by the DR 2 cos 2 8 term in the exponent 
of Eqs. (fTTfl . jT5^ , and (|T3j) . The overlap is minimal at 
these angles as are the resulting amplitudes T^a.- 

The angular momentum content of the scattering ma- 
trix elements T J is obtained by sandwiching the transi- 
tion operators between the initial-state pp and the final- 
state 7r~7r + wave functions, 

T J = / - dRdR'^ 7r (R , )Tfl 2M2 (R,R')*4 =i±1 (R), (27) 



where 0^(R') and ^jE^ (R) still contain the usual an- 
gular dependence in the form of appropriate spherical 
harmonics. Therefore, additional angle-dependent terms 
in 1^2/^2(1^) R-') enhance the contribution of higher an- 
gular momenta J. The resulting increase of contributions 
to higher partial waves J > 1 is the main reason for a bet- 
ter description of the experimental data. We will give nu- 
merical evidence for this in another communication |18| . 
A partial wave analysis of the integral in Eq. (|27|l will 
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show that the richer angular dependence of the R2 tran- 
sition operators considerably increases contributions to 
partial waves J = 2 and higher. 



V. LORENTZ TRANSFORMATION OF 
(ANTI)QUARK SPINORS 

So far, we have solely treated the Gaussian part of 
the pion intrinsic wave functions. An s-wave pion wave 
function in momentum space, omitting flavor and color 
components, is given by 



x u(pi)v(pj) exp 



(28) 



The momenta pi and pj denote the antiquark and quark 
respectively, is the pion momentum, (3 is the param- 
eter that determines the size of the pion and u(pi) and 
v(pj) are the usual quark and antiquark Dirac spinors 
for particles with spin s — ^, 



u(p,m s ) 
v(p,m s ) 



\/E + m 

O-P I Xr, 

) \ I-)" 

V hi + m 



or in terms of helicity spinors 
u(p,A) = 



2p\ Xa(P) 



\ \JE+rn 
/ -2pA 

«(p,A) = >f±^ 
\ VE + ? 



X-x(P)(") 5 ' 



(29a) 
(29b) 

(30a) 
(30b) 



Our approach differs from that of Maruyama et al. 
|13|. who use an approximation of the Is quark wave 
function in the MIT bag in r-space as introduced in 
Ref. 0, 0| . In Ref. the quark wavefunction is 



^(r,) = (l + 3C 2 /2)r* x 



exp{-rf /2i^} ( ) x , 



(31) 



where Rq is the size parameter, which can be related to 
the bag radius. The parameter £ determines the prob- 
ability that the quark is in the lower component of the 
Dirac spinor and is referred to as the small relativistic 
component of the wave function. The value of C in the 
bag model has been estimated to be 0.36 for zero quark 
mass. The authors of Ref. also Lorentz transformed 
the Gaussian part of the wave function in Eq. 131|) into 
the cm. but did not apply this transformation to the 
spinor. In their analysis, Maruyama et al. found that 
the small relativistic components have little influence on 



the predictions of relative ratios of pp annihilation am- 
plitudes into mesons. The effect of including the Lorentz 
contraction, on the other hand, is to enhance the produc- 
tion of two pions. 

In Sec. IIIII the momentum space transition operators 
in Eqs. @, ljTU)l. (|2T|) . and were calculated mak- 
ing use of just the Gaussian part of the wave function 
of Eq. JSSJ. In the present section we investigate how 
Lorentz boosts of quark spinors in the pion wave func- 
tion of Eq. (|28|l affect these calculations. 

In the pion rest frame S w the quark four-momentum is 
p, while in the cm. frame S cnL of the reaction pp — > Tt~ir + 
the quark four- momentum is p' . They are related by the 
inverse Lorentz transformation l^ 1 = ^ _1 (/3) as follows 



P 



p. 



(32) 



Here, (3 — v/c is a vector along the pion momentum P^. 

Helicity states in frame are denoted by |p, A)^ 
while in the cm. frame S cm they are |p',A')s cm . The 
two sets of states are related by 



\p',X)s m = u(r 1 )\p,\)s„ 



(33) 



where U{l~ r ) is a unitary operator effecting the inverse 
Lorentz transformation. It can be shown for any Lorentz 
transformation I that the operation of U(l) on the state 
|p, A) is equivalent to a rotation r(l, p) in spin space. For 
a detailed derivation of Lorentz transformations of helic- 
ity states see for instance Ref. The rotations r(l,p) 
that yield the cm. states |p, A)s cm are the Wick helicity 
rotations. With the appropriate 2s + 1 dimensional ma- 
trix representation T> xx , [r(l, p)], where s = ^ is the spin 
of the (anti)quark, one can rewrite Eq. (|33l) as 



|P,A>, 



\p , A')s m = V tlX r[r(l,p)]\r 1 p,n) 3w . (34) 



A sum over \x is implied. The quark wave functions are 
obtained from 



W(a:)|p,A> 



1 



u(p,A)e~ 



■ IpX 



(35) 



where ip(x) is the quark Dirac field. The correspon- 
dence |p, A) < — > u(p, A) between the physical state and 
its wave function realization yields the desired Lorentz 
transformation of the Dirac spinor 

u(p', A')s cm = Wxx (r)u(r 1 P , A) Sw (36a) 

and similar arguments lead to 

v(p', A')s cm = V s xx , (r-X^p, A) s ^ . (36b) 

For the general case of a boost of the axes I ((3), with 
(3 = (/3, 6p, (j>p), the matrices T> xx , are 



Vi x ,=e"^'^d\ xi {e w ), 



(37) 



where rj is a function of the angles 9p and <pp and the 
Wick angle 8 W 0] is given by 



cos 9 W = — (p + j3E cos 5) . 
V 



(38) 
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Here, E = \/p 2 + m 2 is the energy of one (anti)quark 
as in Eqs. QZOk ) and QZQb ). 7 = E cm /2m n , p = |p|, 
p' = I Z 1 p» I , and < S < it is the angle between p and 
(3, defined by 



cos S — 



P-P = P„ P 

(3p PkP ' 



(39) 



For Lorentz boosts Z z (/3), for which the z-axis is along 
Pjr in the cm. system and Op = (f>p = r) = 0, the expres- 
sion in Eq. H37|l simplifies to 



(40) 



Since s = i, only two matrix elements are required |17| : 



^1/2,1/2 (^ w ) ~ cos 2 



1 1 

— I — cos ( 

2 2 



- + lp(p + 0EcQBg), (41a) 



,1/2 



i/2,i/ 2 ( e w) = sin— = W---cos6» w = 



1 1 



^-^7(p + /3£cos<5). (41b) 



This means that the quark helicity spinors in S cul are 
expressed in terms of Dirac spinors in S„ as 



C\J E + m \ 1 i n 

H-jpl X i (p'K^,v( w) + 



-|r' P | X_ 1 (p')^ / 1 2 /2 ,v( w)(42a) 



and similarly for the antiquark spinors 



< JE+m 

/eT~, 

I'" X P| 



v(p',\') = [ \ x i(p ')^ 2 (0 W ) 

y y/E + m J 5 ' 1 



>§±^- JXxfc'Kv-i/aPw)- (42b) 



To obtain the above expression for antiquark spinors the 
property d\y(—8 w ) — dy\(8 w ) has been used. Note that 
the spinors now depend on Ti a b via the relativistic factor 
7 present in Eqs. I|41al) and (|41b|l as well as in the Lorentz 
boost J -1 . 



VI. RELATIVISTIC SPINOR CORRECTIONS 
TO R2 AND 42 DIAGRAMS 



Having shown in Sec. IIIII how the Lorentz transforma- 
tion of the spatial part of the intrinsic pion wave function 
alters the annihilation operators, we will now discuss the 
effect on the operators due to modifications in the spinors 



of the final quarks i — 1', 2' and antiquarks i = 4', 5' . 
The A2 and R2 spin-matrix elements read schematically 



M s 



{A*}=-§ 



[«2'S4']s„=0 



X V R2 /A2 \[SlS2S 3 ] s 1)® [s 4 s 5 s 6 ] s 1), (43) 



where Vr2/A2 is either one of the four operators intro- 
duced in Eqs. 0, ©, and flSOJl. As before, we 
neglect relativistic corrections to the pp wave functions, 
hence, only the final quark and antiquark helicity states 
are Lorentz transformed. Note that the operators Vr2/A2 
mentioned above are in the cm. frame, where they act 
between initial- and final-state spin wave functions. 

In the derivation of the transition operators in Eqs. J7J, 
(JHJ, i|19l) . an d 1|20|) . only leading order terms in p were 
kept. Here the same strategy is followed: we make use of 
Eqs. I|42a|) and l|42b|) but we do apply these boosts to the 
Pauli spinors in the pions, each consisting of a quark i 
and antiquark j with spin states \sim,i) and \sjirij) (with 
corresponding helicities and Xj) and m, 7^ mj = ±|: 

I^tt = m n = 0) = (^\simi)®\sjmj)-\simj)®\sjrrii 

The Pauli spinors in the cm. frame are expressed in 
terms of the Pauli spinors in S„ by 



X y ( P ') = d\% A , (8 w ) Xl (p') + y (0 W )X t (P') (44a) 

' ' 2 ' ' 2 

for particles and 

Xv(p')-4' 2 i/ 2 ( w)x 1 (p')+4'-i/2^wk AP') (44b) 

• ' 2 1 ' 2 

for antiparticles. 

The Pauli matrices <r in the transition operators 
Vr2/A2 ac t on t ne boosted spinors Xai with i = 
l',2',4',5'. The spin-matrix elements pertinent to the 
^2/A2 operators are thus modified since the rotation ma- 
trix elements d±^ 2 1 ^ 2 (8 W ) in (|41a|l and (|41b(l come into 
play. We make use of the fact that in the cm. frame 
the transverse components of the quark momenta are 
much smaller than the component along the boost di- 
rection (from now on we refer only to the quark but the 
discussion is identical for the antiquark). We can thus 
approximate 



,1/2 



\r\P)p\ = v / ^ 2 + 7 2 b„ +pe) 2 

~ j(p-P„+pE), (45) 

which we write as p' ~ 7(pcos<5 + (3E). This approxi- 
mation is valid except for a small region where cos 6 = 
—f3E/p. With this, the following approximation for the 



matrix elements d 



,1/2 



1/2 

±1/2,1/2^ 



,) holds 



"±1/2,1/2 v 



'1 l p + (3E cos S 
2 2 p cos 5 ■ 



■j3E' 



(46) 
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The numerator and denominator in Eq. (|46() are very sim- 



ilar and the matrix elements d 



1/2 



are functions 



"±1/2,1/2V 

of (3 rather than of 7. 

To illustrate the role of the cos 5 term in Eq. 1(46(1 , we 
consider two cases cos<5 = ±1, i.e. for which the quark 
momentum is parallel or anti-parallel to the boost direc- 
tion (3. Recall that for cos<5 = ±1, of course, p ± vanishes 
in which case Eqs. (|45|l and 146(1 are exact. 

The following relationships exist between helicity and 
momentum orientation of a quark in the cm. frame: 



If cos 5 



d 1/2 (( 
"1/2,1/2^ 

yl/2 / 



= 1 

w) = 



(47) 



In other words, if the quark momentum in the pion rest 
frame is aligned with the boost direction, then helicity is 
conserved under the boost. In this case there is no Wick 
rotation and W = 0. Of course the momentum p' and 
boost vector f3 will also be parallel in the cm. frame. 

However, if momentum in the pion rest frame is anti- 
parallel to the boost direction /3 then helicity in the cm. 
frame can flip, as seen from Eq. 1(48(1 and now 9 W = tt. 



If cos <5 = — 1 



d 1/2 a 

U l/2,l/2^ 

d 1/2 ( 

-1/2, 1/2^ 



,) =0 

w) = 1 



(48) 



The relativistic effect described in Eqs. l|4*7| and l j4"S|l 

can be interpreted as follows. Assume a quark inside 
the pion with a given helicity and its momentum ori- 
ented anti-parallel to the boost direction (3. For very 
large boosts, an observer in the cm. system perceives a 
reversed quark-momentum. As the spin direction is not 
changed, this implies a flipping of the helicity in the cm. 
frame, which confirms that for cos<5 = —1 the matrix 
element d l ^jt, 1 ^ 2 (6' w ) is zero whereas d x J 2 ^ 2 1 i 2 {0 vi ) is one 

as seen in Eq. ((48(1 . On the other hand, if the quark mo- 
mentum in the pion frame is parallel with /3, then helicity 

will not flip, as the quark momentum in the cm. is still 

1/2 

seen in boost direction and d\j 2 1/2 (^ w ) mus t be zero in 
agreement with Eq. 1(471) . Hence, helicities can flip from 
one reference frame to the other depending on the initial 
orientation of the quark momentum with respect to the 
boost direction. 

For intermediate angles —1 < cos<5 < +1, Eqs. I(44a() 
and l(44b() describe the degree of helicity flip due to the 
change of reference frame. Since the rotation matrices are 
proportional to (3, the relativistic effects from the spinors 
tend to be much smaller than those from the spatial com- 
ponents discussed previously in Section III. 



VII. CONCLUSION 



In this paper, the effects of Lorentz transformations on 
intrinsic pion wave functions in the reaction pp — ► n~ir + 



are investigated. This is done within the framework of a 
constituent quark model used previously 0, 0, HI to de- 
scribe da/dSl and Ao n data of this reaction. The Lorentz 
transformations are effected on both the spatial part as 
well as the quark spinor components in the intrinsic pion 
wave functions. We find that the coordinate-space wave 
function is strongly contracted along the boost direction 
and that the spatial part of the annihilation amplitudes 
depends on 7 and j 2 . The spinors are modified by Wick 
rotations and the relativistic effects due to spin are there- 
fore proportional to j3 rather than to the boost factor 7. 

In addition to the non-planar quark rearrangement dia- 
grams R2, we also compute the transition operators orig- 
inating from the planar A2 diagrams for pure 3 S\ as well 
as a mixture of 3 S\ and 3 Pq annihilation mechanisms. 

The non-relativistic annihilation operator for the non- 
planar R2 diagrams was given in Ref. |2j. In this paper 
we obtain new terms in the operator that affect its range 
and introduce additional angular dependence. The de- 
pendence on the cm. energy y/s is now manifest via the 
relativistic boost factor 7. Most of the new relativistic 
factors are considerably larger than the non-relativistic 
ones. Therefore, the geometry of the annihilation ampli- 
tudes is drastically modified. Nonetheless the selection 
rules described in Ref. ?| are preserved. The A2 am- 
plitudes are simpler than the R2 amplitudes and they 
depend only on the relative distance of the proton and 
antiproton. 

We will show [Ts] that the modifications in the an- 
nihilation amplitudes, due to relativistic considerations, 
strongly improve the description of the cross section and 
analyzing power as measured at LEAR p]. Therefore 
with a relativistic treatment one obtains a description 
of the reaction pp — > ir~n + similar or better than with 
the ad hoc increase of the particle radii which was the 
approach in Ref. |7J. 
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